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Abstract: Recent comparisons of minimal supergravity (mSUGRA) model predictions 
with WMAP measurements of the neutralino relic density point to preferred regions of 
model parameter space. We investigate the reach of linear colliders (LC) with yfs = 0.5 
and 1 TeV for SUSY in the framework of the mSUGRA model. We find that LCs can cover 
the entire stau co-annihilation region provided tan/3 ~ 30. In the hyperbolic branch/focus 
point (HB/FP) region of parameter space, specialized cuts are suggested to increase the 
reach in this important "dark matter allowed" area. In the case of the HB/FP region, the 
reach of a LC extends well past the reach of the CERN LHC. We examine a case study 
in the HB/FP region, and show that the MSSM parameters \i and M2 can be sufficiently 
well-measured to demonstrate that one would indeed be in the HB/FP region, where the 
lightest chargino and neutralino have a substantial higgsino component. 
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1. Introduction 



In recent years, supersymmetric models have become increasingly constrained by a variety 
of measurements These include determination of the branching fraction BF(b — ► 
57) [Q], the muon anomalous magnetic moment = {g — 2)^/2 || and most recently, 
the tight restriction on the density of relic dark matter from the Big Bang, as determined 
by the WMAP experiment^]. Analysis of WMAP and other data sets have determined a 
preferred range for the abundance of cold dark matter Q, [| : 

VcDMh 2 = 0.1126lg$8i> 2a level. (1.1) 

Within the minimal supergravity (mSUGRA) framework^], the lightest neutralino 
is usually the lightest SUSY particle. Since i?-parity is assumed to be conserved, this 
neutralino is stable and provides a good candidate for cold dark matter 0. While it is 
possible that the relic density of neutralinos may make up almost all the cosmological 
dark matter, the possibility that dark matter, like visible matter, is made up of several 
components cannot be excluded at this point. In our analysis we will, therefore, interpret 
the WMAP measurement (|1.1| ) as an upper bound, 

h 2 < 0.129, 

on the relic density in neutralinos. The mSUGRA model is characterized by four SUSY 
parameters together with a sign choice, 

mo, rn 1 / 2 , Aq, tan/? and sign(fx). (1-2) 

Here mo is the common mass of all scalar particles at Mqut, mi/2 is the common gaugino 
mass at Mgut, Aq is the common trilinear soft term at Mqut, tan/3 is the ratio of Higgs 
field vacuum expectation values at the scale Mz, and finally the magnitude - but not 
the sign - of the superpotential [i term is determined by the requirement of radiative 
electroweak symmetry breaking (REWSB). In addition, we take mt = 175 GeV. 

The recent measurements of SIcdm^? \ BF(b — ► 57) and have considerably modified 
our expectations for the regions of mSUGRA model parameter space that may be realized 
in nature. Several years ago, relatively low values of mo and mj/j appeared to be preferred 
because these led to sufficiently light sleptons which in turn led to efficient neutralino 
annihilation in the early universe via ^-channel sfermion exchange diagrams; this resulted 
in values of flcDMh 2 < 1, and a universe at least as old as its oldest constituents || ||. 
The low mo and m^ 2 region of parameter space has been dubbed by Ellis et al. as the 
"bulk" annihilation region [p!cf|. The bulk annihilation region was also favored by fine-tuning 



estimates in the mSUGRA model [11]. 



Currently, the bulk annihilation region of the mSUGRA model is disfavored by: 

1. a value of m^ typically below bounds from LEP2, which require m^ > 114.4 GeV in 
the case of a SM-like light Higgs scalar /i[12[; 
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2. a value of BF(b — * 57) that is either above (fi < 0) or below (fi > 0) current 
measurements 1 ; 

3. for larger values of tan j3, the bulk region leads to large negative (fx < 0) or large 
positive (/i > 0) contributions to the muon magnetic dipole moment. Because of 
uncertainty in the SM prediction of an, caution is advised in the interpretation of 
what the E821 data tell us about the existence of physics beyond the SM. 2 

The favored regions of mSUGRA parameter space now include [|14|, [15|, 16 

• the stau co-annihilation region at low mo where m^ x ~ m z 1 > an< ^ where f\ — Z\ and 
fx — fx annihilation in the early universe also serve to reduce the neutralino relic 
density to sufficiently low values [TT 



the ^4-annihilation funnel at large tan/3 where run and tua — 2m^ and Z\Zx 



A, H — > // (/'s are SM fermions) through the very broad A and H resonances! 18], 
and 

the hyperbolic branch/focus point region^, ^] (HB/FP) at large niQ near the 
edge of parameter space where \x becomes small, and the Z\ has a significant higgsino 



component which facilitates a large annihilation rate[|22j, |23|, Bfl. The location of this 
region is very sensitive to the value of ra t . 

The HB/FP region predicts multi-TeV squark and slepton masses so that supersymmetric 
contributions to BF(b — > 57) and a M are suppressed; these measurements are thus expected 
to be close to the SM predictions. Furthermore, the relatively high scalar masses help 
suppress potential flavor changing and CP-violating supersymmetric processes, and offer 
at least a partial decoupling solution to the SUSY flavor and CP problems^]. Finally, 
Feng et al. have shown that the low mx/2 part of the hyperbolic branch has relatively low 



fine-tuning [20 1, in spite of mo being very large. 

Given these expectations for where supersymmetry might lie, it makes sense to re- 
evaluate the prospects for SUSY searches at various collider and dark matter search ex- 
periments. In Ref. |26|| , the reach of the Fermilab Tevatron for isolated trilepton events 
from W\Zi production was extended to very large mo values to include the HB/FP region. 
In the HB/FP region, production cross sections increase due to decreased chargino and 
neutralino masses, but the visible energy from W\ — ► / '/' ' Z\ and Z^ — ► f fZ\ decays also 
decreases, reducing detection efficiency. It was shown that the Tevatron reach in m^/2 for 
isolated 3£ events did in fact increase in the HB/FP region. The reach of the CERN LHC 
was also worked out in Ref. |27j, where it was found that values of mx/2 ~ 1400 GeV 
(~ 700 GeV) could be probed in the stau co-annihilation region (HB/FP region) with an 
integrated luminosity of 100 fb _1 . For higher values of mx/2 i n the HB/FP region, m^ 

A combination of measurements from the ALEPH, BELLE and CLEO experiments yield BF(b —* 
try) = (3.25 ± 0.54) x 10 -4 , while the SM prediction is (3.6 ± 0.3) x 10~ 4 g. 

2 In arecent analysis, Davier et o«.Q find Aa M = (22.1±11.3) x 10~ 10 ((7.4±10.5) x 10~ 10 ) [errors added 
in quadrature] depending on whether the hadronic vacuum polarization is estimated using e + e~ — > hadrons 
(t decay) data. 
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becomes very large while the visible energy from W\ and Z2 decays becomes small, so 
that signal detection becomes difficult. Thus, the high m\/2 part of the hyperbolic branch 
currently seems beyond LHC reach. It is interesting to note that in the HB/FP region of 
the mSUGRA model, the higgsino component of the Z\ is sufficiently large to yield direct 
neutralino-nucleus scattering rates within the reach of Stage 3 direct dark matter detec- 
tion experiments such as Genius, Cryoarray and Zeplin-4|28|. In addition, in the HB/FP 
region one may expect detectable rates for detection of neutrinos arising from neutralino 
annihilation in the core of the sun or the earth, and also large rates for cosmic photons, 



positrons and antiprotons due to DM annihilation in the galactic halo[[23|, 24]. 

In the HB/FP region, since becomes small, charginos are light. This then implies 
that there would be a large rate for chargino pair production at e + e~ linear colliders (LCs) 
operating with center-of-mass energy yfs ~ 0.5 — 1 TeV over most of the HB/FP region. 
However, since the W\ — Z\ mass gap also becomes small, it is not clear that linear collider 
experiments would access the entire kinematically allowed chargino pair production region. 

Our main goal in this paper is to assess the reach of linear colliders for SUSY in the 
mSUGRA model [29, 30 1, paying particular attention to the HB/FP region. We note here 
that previous reach estimates of linear colliders for SUSY in the mSUGRA model extended 
only up to mo values as high as 800 GeV[29|- well below the HB/FP region. We find that 
a linear collider, using standard cuts for chargino pair events, can explore the low m\/2 
portion of the HB/FP region. To explore the high mi/2 part of the hyperbolic branch, new 
specialized cuts are suggested. With these cuts, it appears possible to probe essentially all 
of the HB/FP region (up to my 2 = 1-6 TeV) where charginos satisfy the LEP2 bounds, 
and where chargino pairs are kinematically accessible at a linear collider (LC). This then 
provides the first example of a SUSY parameter space region which is accessible to linear 
e + e~ colliders, while likely remaining out of reach of LHC experiments! This is especially 
interesting since the HB/FP region is one of the three qualitatively different mSUGRA 
parameter space regions allowed by dark matter and other constraints. 

If indeed a SUSY signal from charginos is detected, the next step would be to try 
and determine the associated weak scale parameters: fj,, M2 and tan/3[31|. We explore a 
particular case study in the low m^j region of the HB/FP region, and show that at least in 
this case [i and M2 should be measurable. This measurement would give a firm indication 
of the large higgsino content of the light chargino and Z\ and, together with the fact that 
sfermions are not detected either at the LC or at the LHC, provide a strong indication that 
SUSY in fact lies in the HB/FP region. 

The remainder of this paper is organized as follows. In Sec. 2, we examine the problems 
associated with detecting SUSY in mSUGRA parameter space at an e + e~ LC, with special 
emphasis on the HB/FP region, and show that these can be overcome. We propose novel 
cuts that allow the large m\/2 part of the hyperbolic branch to be explored via chargino 
pair production. We then show our projections for the SUSY reach of a LC operating at 
\fs = 0.5 TeV or \fs = 1 TeV, assuming 100 fb _1 of integrated luminosity. In Sec. 3, we 
compare the reach of a LC with the reach of Fermilab Tevatron luminosity upgrades and of 
the CERN LHC. Furthermore, we also show the dark matter allowed region of mSUGRA 
parameter space to elucidate the performance of these experiments in this preferred region 
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of parameter space. In Sec. 4, we perform a case study in the HB/FP region, where chargino 
pair production is accessible to a y/s = 0.5 TeV LC. The masses and can be 

measured via energy distribution end-points. We find that these measurements, along 
with measurement of the total chargino pair cross section, allow a determination of the 
underlying SUSY parameters \i and M2, although tan/3 remains relatively undetermined. 
In Sec. 5, we present our conclusions. 



2. Reach of a Linear Collider in the mSUGRA model 

In our signal and background computations, we use ISAJET 7.69[[32jl which allows for the 
use of polarized beams, and also allows for convolution of subprocess cross sections with 
electron parton distribution functions (PDFs) arising from both initial state bremsstrahlung 
and also beamstrahlungp3| . We use the ISAJET toy detector CALSIM with calorimetry 
covering the regions —4 < r\ < 4 with cell size A77 x A(p = 0.05 x 0.05. Electromagnetic 
energy resolution is given by AE em / E em = 0.15/ V 'E em © 0.01, while hadronic resolution 
is given by AE^/E^ = 0.5%/^ 0.02, where © denotes addition in quadrature. Jets 
are identified using the ISAJET jet finding algorithm GET JET using a fixed cone size of 
AR = \J Ai] 2 + Ac/) 2 = 0.6, modified to cluster on energy rather than transverse energy. 
Clusters with E > 5 GeV and \r/(jet)\ < 2.5 are labeled as jets. Muons and electrons are 
classified as isolated if they have E > 5 GeV, \rji\ < 2.5, and the visible activity within a 
cone of R = 0.5 about the lepton direction is less than max(Ei/W GeV, 1 GeV). Finally, 
jets originating from 6-quarks are tagged as 6-jets with an efficiency of 50%. 

2.1 Review of previous reach assessment 

The reach of a = 0.5 TeV LC for supersymmetry has previously been evaluated in 



Ref. [29] assuming an integrated luminosity of 20 fb . Reach contours were presented for 
the case of the mSUGRA model in the mo vs. m\/2 plane for Aq = 0, tan (3 = 2 and 10, and 
/i<0. The reach plots in that study were limited to nii/2 < 600 GeV and uiq < 800 GeV, 
i.e. well outside the HB/FP region. 

The region of the mo — mi/2 plane where there should be an observable SUSY signal 
in LC experiments consists of three distinct pieces. 

• At low mo with m 1 / 2 ~ 300 — 500 GeV, slepton pair production occurs at large 
rates. The signal is a pair of opposite sign/same flavor leptons plus missing energy. 



Tsukamoto et al.^Vj suggested cuts of i) 5 GeV < E(£) < 200 GeV, it) 20 GeV 
< E vis . < yfi- 100 GeV, in) \m{ll) - M z \ > 10 GeV, iv) | cos 0(^)1 < 0.9, v) 
—Qicos0£ < 0.75, vi) 9 aC op > 30°, vii) E™ ls > 25 GeV and viii) veto events with 
any jet activity. Here, 6 acop = 7r — cos +pypy)- The reach was evaluated 

by running with right-polarized electron beams where Pi{e~) = —0.9. The beam 
polarization maximizes £r£r pair production, while minimizing background from 
W + W~ production. 

At low mi/2 values, chargino pair production occurs at a large rate. To search for 
chargino pairs, one may look for 1£ + 2j + E mts events. Following Ref. [31], it 
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was required in Ref. [29] to have one isolated lepton plus two jets with i) 20 GeV 
< E ms < v/£ - 100 GeV, it) if E j5 > 200 GeV, then m(jj) < 68 GeV, Hi) E™ is > 
25 GeV, iv) \m{lv) - M w \ > 10 GeV for a W pair hypothesis, v) \cos6{j)\ < 0.9, 
| coe 0(4) | < 0.9, -Q £ cos0 £ < 0.75 and Qicos9(jj) < 0.75, m) aC pp(WW) > 30° for 
a pair hypothesis. The reach for l£+2j+ Et events from chargino pair production 
was evaluated using a left polarized beam with Pl = +0.9. 

• Finally, there exists a small region around mo ~ 200 — 500 GeV and m\/2 ~ 300 — 
350 GeV where neither slepton pairs nor chargino pairs are kinematically accessible, 
but where e + e~ — > Z1Z2 is. In this case, the decay Z2 — * Z\h was usually found 
to be dominant. Since h —* bb with a large branching fraction, bb+ Et events 
were searched for with two tagged b-jets, Ej! %s > 25 GeV and 30° < A^ 6 g < 150°. 
Imposing a missing mass cut tjh > 340 GeV eliminated almost all SM backgrounds, 
so that a signal cross section of 10 fb would yield 10 signal events for integrated 
luminosity 20 fb -1 , where the signal efficiency was found to be 6%. Beam polarization 
of Pl = +0.9 was used. 

In the following, we refer to these as the "standard cuts" . 

The ultimate reach contours found in Ref. [29] generally track the boundary of the 



kinematically allowed regions for W\ and cr pair production. An exception occurs at low 
mo values where selectron pair production is dominant, but where mg R ~ rn Zi m Then, 
the mass gap mg R — becomes so small that there was very little visible energy in 
the slepton pair events, resulting in very low detection efficiency, causing a turnover in the 
reach contours. For tan/3 values larger than those explored in Ref. f29|| , the large r Yukawa 
coupling makes fi significantly lighter than e.R so that close to the boundary of excluded 
region at small thq e^e/j events may still be observable, while f\T\ events are not. 

2.2 Updated reach results 

In this section, we update our earlier reach projections [p9f| for a linear collider. We present 
reach projections for linear colliders with yfs = 0.5 and 1 TeV of energy in the CM frame. 
We also expand the range of mi/ 2 (to 1.6 TeV) and mo (to 8 TeV) beyond the values 



presented in Ref. [29]. This allows us to explore the entire stau co-annihilation strip, the 
^-annihilation funnel and the HB/FP region. In our analysis, we restrict ourselves to the 
trilinear SSB term Aq = 0. For the most part, our results are qualitatively insensitive to 
variations in Aq. An exception occurs for particular Aq choices which may greatly reduce 
the value of m^, and lead to a top squark-neutralino co-annihilation regional. 

Our first results are presented in Fig. |l|, where we show the mo vs. my 2 plane for 
tan/3 = 30, \i > and Aq = 0. The left-most red region at low mo is disallowed because 
the LSP would become a stau. The right-most red region (large mo) is mainly excluded 
by a lack of appropriate REWSB, although this includes as well points with no convergent 
RGE solution as generated by IS A JET. The precise location of the boundary of the large 
mo red region depends somewhat on the computer code, and also on the assumed fermion 
masses. The lower yellow region is excluded by LEP2 chargino searches, which require 
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> 103.5 GeV. In addition, the region below the yellow contour gives a light (SM-like) 
Higgs boson with mh < 114.4 GeV. 

Using the standard dilepton cuts as described above, we find a SM background level 
of a SM = 1-79 fb (0.045 fb) for a LC with yfe = 0.5 TeV (1.0 TeV) with right polarized 
beams using P L (e~) = -0.9. 3 It has been shown in Ref. |3l]] that backgrounds from 2 — * 3 
processes such as e + e~ — > uuZ, e + e~Z or e^vW^ production or 2 — ► 4 processes such as 
e + e~ — » W + W~ production are also efficiently removed by these cuts, at least for yfs = 
500 GeV. Although the cross sections for these processes grow with energy, we expect that 
these cuts will remove the bulk of these backgrounds also at yfs = 1000 GeV; for instance, 
the dominant portion of the eeWW background that comes from "77" collisions will be 
removed by the acoplanarity and other cuts; much of the remaining cross section will have 
p T (WW) < My/ and will also be reduced, though not eliminated by these cuts. However, 
since these backgrounds have not been included in our evaluation it is possible that the 
statistical significance of the signal may be somewhat over-estimated for yfs = 1000 GeV. 
The regions to the left of the lower (upper) blue contour yield a supersymmetric signal at 
the 5(j level assuming 100 fb" 1 of integrated luminosity at the 0.5 (1) TeV LC. An increased 
reach for slepton pairs may be obtained by searching for ditau events originating from stau 
pair production. While we did not perform a detailed signal-to-background analysis for 
this channel, we do show the kinematic limit for stau pair production by a dashed light- 
blue contour, which for tan/3 = 30, lies somewhat above the blue dilepton reach contour. 
Thus, a significant gain in reach might be acquired by searching for acollinear di-tau events 
arising from stau pair production. In particular, this signal may give access to the large 
m i/2 part of the stau co-annihilation region. 

The green contour denotes the reach of a 0.5 or 1 TeV LC for SUSY via the l£+2j+ $t 
channel arising from chargino pair production, using the standard cuts given above. Unlike 
Ref.||], we use P L = 0, and find a background level of 15.5 fb (2.1 fb) for ^/i = 0.5 (1) TeV 
(the beam polarization is not important for this reach contour). For most of parameter 
space, the reach contours follow closely along the = 250 (500) GeV mass contours, 
indicating that chargino pair production can be seen with standard cuts almost to the 
kinematical limit for chargino pair production. An exception occurs when mo becomes 
very large, in the HB/FP region. Around tuq ~ 4000 GeV (mo ~ 6000 GeV) for y/s = 0.5 
TeV (1 TeV), the reach contour departs from the kinematic limit. The termination of the 
reach contour occurs because in this region, the superpotential parameter fx becomes very 
small, and the light chargino W\ and neutralino Z\ become higgsino-like, and increasingly 
mass degenerate. The Q-value from W\ — ► Z\f f decay (the fs are light SM fermions) 
becomes very small, and very little visible energy is released by the chargino decays. This 
causes the detection efficiency for 11 + 2j+ fir events to decrease sharply, leading to a 
corresponding reduction in the reach using the standard cuts. 

To understand what is happening in the HB/FP region, we show relevant sparticle 
masses in Fig. [2|a) for mi/ 2 = 225 GeV, tan/3 = 30, Aq = and \i > versus the 
parameter mo- As mo varies from 1400 GeV to nearly 2200 GeV, i.e. as we approach 

3 In our assessment of SM backgrounds, we have evaluated only the backgrounds from 2^2 processes. 
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Figure 1: Reach of a linear collider for supersymmetry in the mSUGRA model for ^/s = 500 and 
1000 GeV, for tan/3 = 30, A — and \i > 0. The reach via slepton pair production is denoted by 
the blue contour, while standard cuts for chargino pair production yield the green contour. Special 
chargino pair cuts yield the black contour in the HB/FP region. The red region is theoretically 
excluded, while the yellow region is excluded by LEP2 measurements. Below the yellow contour, 
m h < 114.4 GeV. 



the HB/FP region for fixed mi/2 with increasing mr,. As mo increases, \n\ is seen to be 
decreasing. Since the value of SU(2) gaugino mass M2 is essentially fixed, the various 
chargino and neutralino masses also decrease, with the lighter ones becoming increasingly 
higgsino-like. The plot is terminated when the LEP2 limit Tn^ > 103.5 GeV is reached. 
Of great importance is that the W\ — Z\ mass gap is also decreasing, although in this case 
it remains substantial out to the edge of parameter space. 

In Fig. |2|6) , we show the total cross section for various chargino and neutralino produc- 
tion reactions versus mo as in frame a). At the lower mo values, a(W^W^~) pair production 
dominates the total SUSY production cross section. As mo increases, and decreases, 
the other charginos and neutralinos become light as well, and many more reactions "turn 
on" in the HB/FP region. Although we will focus mainly on W^W^ pair production, it is 
important to note that many SUSY production reactions can occur in the HB/FP region, 
and can lead to an assortment of SUSY events from the production and cascade decays of 
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the heavier chargino and neutralino states. 
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Figure 2: Plot of a), sparticle masses and b). sparticle pair production cross sections versus mo 
in the HB/FP region for m 1/2 = 225 GeV, tan j3 = 30, A Q = and ^ > for a = 500 GeV e+e" 
collider. 

A similar plot is shown in Fig. ^ except this time for m^2 = 900 GeV, i.e. in the 
upper regions of the hyperbolic branch. In this case, M<i is much larger than the case 
shown in Fig. ^, and so the heavier charginos and neutralinos remain inaccessible to a LC. 
As mo increases, again and decrease. 4 But in this case the W\ — Z\ mass gap is 
much smaller, reaching only several GeV at the limit of parameter space. Clearly, in this 

4 Since \^\ decreases very rapidly as mo increases and approaches the theoretical boundary of the HB/FP 
region, its precise value is not easy to obtain using numerical methods. The value of fj,, of course, directly 
affects the chargino and neutralino masses. In this figure, we have smoothed out what appeared to be 
rather large numerical fluctuations in two of the bins. 
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4680 4700 4720 4740 4760 4780 4800 4820 

m (GeV) 

b) mSugra with m 1;2 = 900 GeV, tanp = 30, A = 0, n > 



o(all) 



- (W,W,) 



10'- 



a(Z,Z 2 ) 



4680 4700 4720 4740 4760 4780 4800 4820 4840 4860 

m (GeV) 



Figure 3: Plot of a), sparticle masses and b). sparticle pair production cross sections versus mo 
in the HB/FP region for m 1/2 = 900 GeV, tan (3 = 30, A = and fi > for a <fs = 500 GeV e+e~ 
collider. 



upper mi/2 region of the hyperbolic branch, there will be little visible energy emerging from 
chargino 3-body decays, making detection of chargino pair events difficult using standard 
cuts. In addition, as shown in frame b), only W^W^ and Z\Z2 pair production occur, 
so fewer anomalous events are expected in the upper HB/FP region. Since ~ m Wi'' 
there will also be little visible energy from Z2 — > Z\ff decay, so that Z1Z2 production 
will also be more difficult to observe. In the deep HB/FP region where \fj\ <C Mi^, one of 
the neutralinos is mainly higgsino-like with roughly equal components of h u and h^. The 
other neutralinos, being orthogonal to these, thus either have equal magnitudes for their h u 
and hd content, or this content is small. In either case, the ZZ^Z^ coupling is dynamically 
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suppressed |5]J in this region. This accounts for the strong suppression of a(Z 2 Z2) (recall 
that the electron sneutrino is very heavy) in Fig. ^b. 

Coincidentally, the reach of a y/s = 0.5 TeV LC in the U + 2j+ Et channel using the 
standard cuts terminates in the HB/FP region at nearly the same m^ 2 value as does the 



reach of the CERN LHC shown in Ref. [27], and again in Sec. 4 of this paper. Meanwhile, 
the contour of chargino pair kinematic accessibility extends to much higher mi/2 values, 
along the hyperbolic branch. This motivated us to examine strategies to extend the reach 
of a LC to the large mi/2 part of the HB/FP region. 

To find better suited signal selection cuts for the HB/FP region, we examine a par- 
ticular case study for the mSUGRA point which is beyond the projected reach of the 



LHC|7|: 



m , m 1/2 , A , tan/3, sign(n) = 4625 GeV, 885 GeV,0, 30, +1, (2.1) 

for which various sparticle masses and parameters are listed in Table [l|. We will refer to 
this as case 1. Not only is this point inaccessible at the LHC, but most of the sparticles are 
also inaccessible to a LC, with the exception being the lighter charginos and neutralinos. 
While niyy = 195.8 GeV, so that W^W^ pair production occurs at a large rate at a 
yfs = 0.5 TeV e + e~ collider, the W\ — Z\ mass gap is only 14.2 GeV, so little visible 
energy is released in chargino pair production events. 

With this in mind, we generate SUSY events for this case study using ISAJET 7.69 
for a linear collider with ^/s = 0.5 TeV and unpolarized beams, including bremsstrahlung 
and beamstrahlung for background events. The beamstrahlung parameters, defined in Ref. 
|33|] , are taken to be T = 0.1072 with beam length a z = 0.12 mm. 

In Fig. |i]a) we show the distribution of E v i s n,i e from events with 1-lepton and two 
jets expected at a -y/i = 500 GeV LC. The solid black histogram represents the SUSY 
case study, which peaks at very low -EWsiMej as expected from the low energy release 
from chargino decays. The small number of events around E V i S n,i e = 250 GeV is from 
Zh production. The green histogram shows the sum of all 2 — » 2 SM backgrounds 
as generated by ISAJET. The large E v i s ibi e component of these arises from WW, ZZ 
and ti production, where some energy is lost due to associated neutrino emissions. The 
SM background distribution extends to low E v i s m e values, and has a visible shoulder at 
E V isMe ~ due to processes such as e + e~ — > Z — > bb, cc, where the Z can be made by 
convoluting the subprocess reaction with the electron PDF, and the lepton from the decay 
of the heavy flavor is accidently isolated. The bulk of the 2^2 SM background can be 
eliminated by requiring low values of E V i S M e . In this case we require 

20 GeV < E visible < 100 GeV. (2.2) 

The upper limit is chosen to be well above the case 1 signal distribution endpoint to 
accommodate later scans over all mSUGRA parameter space, including points which allow 
larger W\ — Z\ mass gaps, and somewhat harder Edible distributions. We also show a 
red histogram which shows the results of the evaluation of the background from e + e~ — > 
e + e~cc, e + e~bb processes when both initial leptons escape detection when being scattered 
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Figure 4: Distribution in a) E^aMe f° r mSUGRA signal (black histogram) with 
(mo, m 1 /2, Aq, tan/3, sign(^) — 4625 GeV, 885 GeV, 0, 30, 1) after cuts in the first row of Ta- 
ble^. We take \fs = 500 GeV and an integrated luminosity of 100 fb -1 , and adopt beamstrahlung 
parameters T = 0.1072 and o z — 0.12 mm. The ISAJET SM background is shown by the green 
histogram, while the background from 77 — > cc, bb is shown in red. In 6), we show the distribution 
in transverse plane dijet opening angle requiring, in addition, that 20 GeV< f^nstUe < 100 GeV. 
In c), we show the distribution in m(£j), after the aditional requirement cos 4>(jj) > —0.6. The jet 
entering the m(£j) distribution is the one that is closest in angle to the lepton direction. 



at a very small angle. We evaluated this background, which mainly arises from photon 
photon collisions, using the PYTHIA event generator |36[| . In the case of bb production, the 
isolated lepton arises from semi-leptonic b — > civ decay, and the jets come one from a b 
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quark, and the other from the charm quark. 

In 77 — ► bb events, the b and b will typically emerge 
back-to-back in the transverse plane. Thus, in Fig. ||&) 
we plot the distribution in cos <p(jj), where 4>(jj) is the 
transverse dijet opening angle. The signal is distributed 
over a range of cos (f>(jj) values, and actually peaks at 
cos (f>(jj) ~ 1. The background peaks at cos <j>(jj) ~ 
— 1, so we require a cut of 



cos (f>(Jj) > —0.6. 



(2.3) 



Finally, any surviving background arising from bb 
or cc production followed by semileptonic heavy fla- 
vor decay is likely to have a jet-lepton invariant mass 
bounded by the heavy flavor mass (at least at parton 
level). In Fig. ||c) we show the distribution in m(£j) 
where we form the invariant mass from the jet which 
is closest to the isolated lepton in space angle. Some 
additional background removal at low cost to signal is 
gained by requiring 



parameter 


value (GeV) 


M 2 


705.8 


Mi 


372.2 


/' 


185.9 


m-g 


2182.7 




4893.9 


m iL 


4656.1 


Wi 


195.8 


m w 2 


743.5 


m z 1 


181.6 




196.2 


^3 


377.3 


m z A 


760.0 


m A 


3998.3 


m h 


122.0 




0.0104 


BF(b -> 87) 


3.34 x 10~ 4 




0.6 x 10~ 10 



m(£j ne 



> 5 GeV. 



(2.4) 



At this point, the distribution is clearly dominated 
by signal. The cross sections in fb after each cut for 
signal and background are shown in Table |2], where we 
include in addition background from the 2 — > 4 process 
£uqq' (evaluated using CompHEP ]37fl) . In the 



Table 1: Masses and parame- 
ters in GeV units for case 1 for 
m o, fi/2, Ah tan/3, signal = 
4625 GeV, 885 GeV, 0, 30, +1 in the 
mSUGRA model. The spectrum is 
obtained using ISAJET v7.69. 



e 1 e 

2^4 calculation we eliminate Feynman diagrams such 

as WW pair production which are already accounted for as 2 — > 2 processes in ISAJET. 
The 2 — ► 4 processes are negligible after the E v i s ibie cu t- For all frames of Fig. || we have 
assumed 100 fb -1 total integrated luminosity. We recognize that the dominant background 
to the rather soft signal which comes from two photon collisions will be sensitive to the 
beamstrahlung parameters, and hence to the shape of the beams. Nevertheless, the large 
signal to background ratio that we find encourages us to expect that our conclusions about 
the large reach of LC will be qualitatively unaltered. 

We now require a 5a signal for SUSY events above the total SM background as listed 
in Table ^, for 100 fb -1 of integrated luminosity and scan mSUGRA points in the HB/FP 
region of Fig. |l|. The new result is the black contour, below which the mSUGRA parameter 
space is accessible by the LC at 5a level. One can see that this contour pushes the reach of 
the LC to much higher values of min- The contour peters out at low mi/j values, where 
the visible energy arising from chargino pair production is typically much higher than the 
100 GeV maximum required by our cuts. However, this low m^j region is already well 
covered by the standard chargino search cuts listed at the beginning of this section. The 
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cuts 


case 1 


ISAJET BG 


77 — > cc, bb 


tvqq' 


7], E, AR 


16.2 


897.1 (483) 


9.2 (6.2) 


448 (712) 


20 GeV < E vis < 100 GeV 


14.4 


12.6 (3.5) 


5.4 (4.9) 


0.16 (0.08) 


cos (f>(jj) > —0.6 


13.5 


0.34 (0.2) 


1.1 (1.1) 


0.04 (0.02) 


m(£j) > 5 GeV 


12.9 


0.17 (0.1) 


0.8 (0.8) 


0.04 (0.02) 



Table 2: Cross section after cuts in fb for mSUGRA case 1 signal and ISAJET SM backgrounds, 
two photon background 77 — > cc, bb and the 2 — ► 4 process e + e~ — » Iveqq'. We take y 7 ! = 0.5 TeV 
collider CM energy. The corresponding background for s/s = 1 TeV case is listed in parenthesis. 

new cuts for the far HB/FP region work for e + e _ colliders with \fs = 1 TeV as well. The 
y/s = 1 TeV reach contour extends even beyond the limits of parameter space shown in 
Fig. I 

To complete our SUSY reach contours, we also examined the reach of a LC for SUSY 
via the e + e~ — > Z\Z 2 reaction in the non-HB/FP part of parameter space, where neither 
chargino pair production nor slepton pair production is kinematically accessible. In this 
case, we follow Ref. p9fi in asking for bb+ Et events from e + e~ — ► Z\Z 2 production 
followed by Z 2 — > Z\h, where h — > bb. In Ref. [p9f| , no background was found after a series 
of cuts listed at the beginning of this section. Here, we do not perform complete event 
generation at every point in parameter space, but instead require that 

a(e + e~ -> Z X Z 2 ) x BF(Z 2 -> Z x h) > 2 fb, (2.5) 

which should yield ~ 10 signal events for 100 fb -1 of integrated luminosity assuming an 
efficiency of 5 — 6% as found in Ref. [^]. The resulting reach contour is shown in Fig. [I] 
as the black contour linking the slepton pair reach to the chargino reach contour. It gives 
some additional parameter space reach to a LC, although it is in a dark matter disfavored 
region of parameter space (unless tan /3 is large, and the H, ^4-annihilation funnel cuts 
through it). There is a turnover in the e + e~ — ► Z\Z 2 reach contour at low mo; this occurs 
because in this region, Z 2 — > t\t decay becomes accessible, resulting in a suppression of 
the Z 2 — ► Z\h branching fraction. However, the Z 2 — > f\f signal may also be detectable; 
if so, the gap caused by the turn-over just mentioned would be filled. 

Finally, we show the kinematic limit for e + e~ — * ZH (green-dashed or dotted con- 
tours), e + e~ — ► Ah (orange dashed or dotted contours) and e + e~ — > H + H~ (purple 
dashed or dotted contours). The dashes are for \/s = 0.5 TeV, while dotted are for \fs = 1 
TeV. For tan (3 = 30, these contours always lie below the sparticle reach contours, so if 
heavier SUSY Higgs bosons are seen (at least in the channels mentioned above), sparticles 
should also be seen if SUSY is realized as in the mSUGRA framework. 

In Fig. [| we show LC reach contours for the same mSUGRA parameter plane as in 
Fig. except this time tan/3 = 10. Qualitatively, many of the reach contours are similar 
to the tan/3 = 30 case of Fig. |l|, and in particular, the new cuts designed to access the 
far HB/FP region again allow the LC reach to extend into the high mi/2 section of the 
hyperbolic branch. One difference for tan j3 = 10 results is that the f^f^ kinematic reach 
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Figure 5: Reach of a linear collider for supersymmetry in the mSUGRA model for y/s = 500 and 
1000 GeV, for tan — 10, Aq = and /i > 0. The colors on the various regions and on the different 
contours are as in Fig. [|. 



contour now lies nearly atop the selectron/smuon reach contour using the dilepton cuts 
described at the beginning of this section. 

In Fig. ^, we show the same reach contours in the mo vs. mi/2 plane, but this time 
for tan/? = 45 and \i < 0. The standard slepton pair and chargino pair production reach 
contours are similar to the low tan (3 cases. In this case, the far HB/FP region cuts allow the 
reach to be extended, although the reach contour terminates in the part of the red region 
where the numerical solutions to the renormalization group equations do not converge as 
per the criteria in IS A JET. A tiny region remains inaccessible to our new cuts, between the 
black solid contours and the dashed contours that depicts the kinematic limit for chargino 
pair production. Another feature of the tan/3 = 45 plot is that the stau pair kinematic 
region has expanded even more beyond where signals for cr and /2r may be accessible in 
the dilepton channel. In addition, for this large value of tan/3, the heavy Higgs bosons are 
much lighter than the low tan/3 cases and now there exist regions of parameter space 
where ZH and Ah production may be accessible, while sparticles are not: indeed such a 
situation would point to a large value of tan /3 which would of course be independently 
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Figure 6: Reach of a linear collider for supersymmetry in the mSUGRA model for y/s = 500 and 
1000 GeV, for tan/3 = 45, A = and fi < 0. The colors on the various regions and on the different 
contours are as in Fig. [|. 



measureable from the properties of the detected Higgs bosons [^9|. These regions all occur 
well below the reach of the LHC for SUSY, which will be shown in Sec. 4. Thus, if nature 
chooses tan (3 = 45 in an mSUGRA-like model, and the LC sees only Higgs bosons beyond 
the SM, it is likely that the existence of SUSY would already have been established by 
LHC experiments. 

In Fig. [7], we show our final LC reach plot, taking tan f3 = 52 with [i > 0. In this case, 
there is still substantial reach for chargino pairs via the standard cuts, and the special 
cuts for the far HB/FP region again allow extended reach into this area. For this large a 
value of tan /?, the dilepton reach contour has been completely consumed by the expanding 
forbidden region on the left where f\ becomes the LSP. In addition, the e + e~ — * Z1Z2 reach 
region has shrunk due to the increased branching fraction for Z2 — * fyr decay. However, 
the kinematic reach for stau pairs has greatly increased, and becomes especially important 
for very large tan/?, especially if nature has chosen to reduce LSP dark matter from the 
early universe via co- annihilation with staus. 
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Figure 7: Reach of a linear collider for supersymmetry in the mSUGRA model for y/s = 500 and 
1000 GeV, for tan/3 = 52, A — and [i > 0. The colors on the various regions and on the different 
contours are as in Fig. [|. 



3. Comparison of LC reach with Tevatron, LHC and fi^ /i 2 

In this section, we present an overview of the reach of a LC in comparison to the reach for 
sparticles that can be obtained by the Fermilab Tevatron and the CERN LHC. In addition, 
we show regions of relic neutralino dark matter density in accord with the recent WMAP 
measurements. 

Our first results are shown in Fig. || which shows the same parameter space plane as in 
Fig. |5| In this case, however, we plot the composite reach plot of a yfs = 0.5 and 1 TeV LC 
for discovery of sparticles assuming 100 fb _1 of integrated luminosity. The LC reach plots 
now consist of combined 1.) slepton pair reach via dileptons (we also show the kinematically 
accesible stau pair production contour, if it is beyond the slepton reach), 2.) the chargino 
pair reach via \l + 2j'+ ft? events, with either the standard cuts or the cuts specialized 
for searches in the HB/FP region, and 3.) the region of Z1Z2 — > bb+ $t- In addition, we 
superimpose on this plot the reach of the Fermilab Tevatron for SUSY via the clean trilepton 
signal originating from pp — > W1Z2X — > 31+ $t + X, where X denotes assorted hadronic 
debris. The Tevatron reach was extended into the HB/FP region in Ref. [26]; we show the 
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optimistic reach assuming a 3a signal with 25 fb" 1 of integrated luminosity. In addition, we 
show the reach of the CERN LHC as derived in Ref. [27|, assuming 100 fb _1 of integrated 
luminosity. We have also added to the plot the green region, which denotes parameter space 
points with relic density ^^Ji? < 0.129, as required by the recent WMAP measurements. 5 



For our relic density calculation pO|] we have evaluated all relevant neutralino annihilation 
and co-annihilation processes in the early universe using the CompHEP program. We have 
implemented relativistic thermal averaging of the annihilation cross section times velocity, 
which is useful to get the appropriate relic density in the vicinity of s-channel poles (where 
the annihilating neutralinos may have substantial velocities) using the formulae of Gondolo 
and Edsjof4lJ. 

The dark matter allowed region splits into three distinct regions for tan/3 = 10. On 
the far left of the plot at low tuq is the stau co-annihilation region, which blends into the 
bulk annihilation region at low values. Note that the bulk region is largely below the 
LEP2 rrih = 114.4 GeV contour. We also see that the stau co-annihilation region extends 
to mi/2 values as high as ~ 900 GeV. For tan/3 = 10, we see that a = 0.5 TeV e + e~ 
collider should be able to scan much of the stau co-annihilation region, while a y/s = 1 
TeV machine can cover it entirely (as can the LHC). Another region of relic density is the 
small strip at constant mi/2 ~ 100 GeV, where neutralinos can annihilate through the 
narrow s-channel pole from the light Higgs boson h. This region can be covered by all 
the colliders, including the Fermilab Tevatron. Finally, adjacent to the REWSB excluded 
region at large mo is shown the dark matter allowed region in the HB/FP region, where the 
LSP has a significant higgsino component, which facilitates neutralino annihilation to WW 
and ZZ pairs in the early universe. The Fermilab Tevatron reach does not extend into 
this regime. The ^fs = 0.5 TeV LC can explore the kinematically allowed portion of the 
lower HB/FP region (which is the region favored by the fine-tuning analysis of Ref. |pCf| ) 
via standard cuts and the cuts specialized to the far HB/FP region. The yfs = 1 TeV 
colider can explore all the HB/FP region which is dark matter allowed, until rnwj becomes 
greater than ~ 900 GeV. The portion of the HB/FP region with m^ 2 > 900 GeV, while 
allowed by dark matter as well as other experimental constraints, becomes more difficult 
to reconcile with fine-tuning considerations. The new cuts proposed in Sec. 2 allow the 
LC SUSY search region to extend well beyond the reach of the CERN LHC, which extends 
only to mi/2 ~ 700 GeV. The LHC reach is limited in the high m^ 2 part of the hyperbolic 
branch because sfermions and gluinos are too heavy to be produced at an appreciable 
rate. Chargino and neutralino pairs can still be produced at the LHC in the high mi/2 
part of the hyperbolic branch, but the soft visible energy emanating from chargino and 
neutralino decay makes detection above background very difficult. The high m\n part of 
the hyperbolic branch yields a first example of a region of mSUGRA model parameter space 
where sparticles can be discovered at a LC, whereas the CERN LHC reach for sparticles has 
petered out. Moreover, this additional reach area comes in precisely at a very compelling 
dark matter allowed region of the mSUGRA model. 

In Fig. [9|, we show the same plot, except this time for tan/3 = 30. Many features of the 

The WMAP allowed region including the lower bound would appear as a very narrow strip following 
the border of the green region. 



-17- 



with tanfi = 10, A Q = 0, jlx > 
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Figure 8: Reach of a y/s = 0.5 and 1 TeV LC for sparticles in the mSUGRA model for tan/3 = 10, 
Aq = and n > 0. We also show the reach of the Fermilab Tevatron assuming 10 fb _1 of integrated 
luminosity (for isolated trileptons) and the reach of the CERN LHC assuming 100 fb _1 of integrated 
luminosity. Finally, the green shaded region shows points where the relic density /i 2 < 0.129 as 
dictated by WMAP. 



plot are qualitatively similar to the tan (3 = 10 case. In this case, the stau co-annihilation 
corridor now extends up to m\/2 values as high as 1050 GeV. The entire stau co-annihilation 
corridor can potentially be explored by a yfs = 1 TeV LC, but only if a stau pair search is 
made in addition to the dilepton search. In this case, the Tevatron reach extends just to 
the tip of the dark matter allowed HB/FP region. The LHC reach in the HB/FP region is 
again limited to mx/2 < 700 GeV values, while the \fs = 0.5 and especially the y/s = 1 TeV 
LC can explore much of the HB/FP region, even for mi/2 values far in excess of 700 GeV. 

The dark matter relic density is qualitatively different for the case of tan j3 = 45, /i < 
shown in Fig. |lfj. Here, a large new dark matter allowed region has emerged, namely the 
yl-annihilation funnel which is characteristic of the mSUGRA model at very large tan (3. 
As tan/3 increases, the derived value of txia decreases, until a region where txia — ^ m z 
arises, where neutralinos can efficiently annihilate through the very broad A and also the H 
s-channel resonances. It can be seen from the figure that the ^-annihilation funnel region 
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Figure 9: Reach of a y/s = 0.5 and 1 TeV LC for sparticles in the mSUGRA model for tan (3 = 30, 
Aq = and n > 0. We also show the reach of the Fermilab Tevatron assuming 10 fb _1 of integrated 
luminosity (for isolated trileptons) and the reach of the CERN LHC assuming 100 fb _1 of integrated 
luminosity. Finally, the green shaded region shows points where the relic density ^j^ 2 < 0.129 as 
dictated by WMAP. We denote the kinematic limit for stau pair production at LCs by a dashed 
black contour. 



extends well beyond the reach of both the \/s = 0.5 and 1 TeV LC. In addition, the stau 
co-annihilation strip rises to m\/2 values that are also beyond the reach of a 1 TeV LC. 
The CERN LHC can explore essentially all of the ^-annihilation funnel for this particular 
value of tan/3 and sign of /x. Also, in this case, the yfs = 0.5 TeV LC can explore the 
HB/FP region only up to m^j ~ 600 GeV where the = 250 GeV contour intersects 
the excluded region. The 1 TeV LC has a reach that extends again well beyond the limit 
of the LHC reach in the HB/FP region. 



In Fig. 11, we show the tan (3 = 52 mSUGRA plane for \i > 0. In this case, the effect of 
the yl-annihilation funnel is just beginning to enter the mo vs. mi/2 plane from the left, so 
that points along the low rriQ forbidden region have a low relic density because neutralinos 
can annihilate via stau coannihilation, via t-channel slepton (mainly stau) exchange (low 
m i/2) an d partly due to annihilation through the s-channel A resonance. In this case, the 
A resonance corridor is actually off the plot, but since the A width is so large (r^ ~ 25 GeV 
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with tan$ = 45, A Q = 0, ju, < 
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Figure 10: Reach of a y/s — 0.5 and 1 TeV LC for sparticles in the mSUGRA model for tan (3 = 45, 
Ao = and /z < 0. We also show the reach of the Fermilab Tevatron assuming 10 fb _1 of integrated 
luminosity (for isolated trileptons) and the reach of the CERN LHC assuming 100 fb _1 of integrated 
luminosity. Finally, the green shaded region shows points where the relic density ^^ 1 h 2 < 0.129 as 
dictated by WMAP. We denote the kinematic limit for stau pair production at LCs by a dashed 
black contour. 



for mi/2 ~ 600 GeV), the value of 2m^ can be a few partial widths away from resonance 
and still give significant contributions to the neutralino annihilation rate. For this large a 
tan/3 value, the stau co-annihilation strip reaches mi/2 values far beyond the reach of LCs 
or even the LHC. In the HB/FP region, the new cuts presented in Sec. 2 again give the 
LCs a reach well beyond the LHC for rrii/ 2 > 700 GeV. 

4. Determination of Model Parameters in the HB/FP Region 

Once a signal for supersymmetry is established at a LC, then the next task will be to 
scrutinize the signal to elucidate production and decay processes, extract sparticle masses, 
spins and other quantum numbers, and ultimately to determine parameters of the under- 



lying model. Many groups have examined different case studies [42]. In this section, we will 
examine a case study in the low m^ 2 part of the HB/FP region in an attempt to extract 
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Figure 11: Reach of a y/s — 0.5 and 1 TeV LC for sparticles in the mSUGRA model for tan (3 = 52, 
Ao = and n > 0. We also show the reach of the Fermilab Tevatron assuming 10 fb _1 of integrated 
luminosity (for isolated trileptons) and the reach of the CERN LHC assuming 100 fb _1 of integrated 
luminosity. Finally, the green shaded region shows points where the relic density ^^ 1 h 2 < 0.129 as 
dictated by WMAP. We denote the kinematic limit for stau pair production at LCs by a dashed 
black contour. 

the underlying parameters of the MSSM, which may in turn point to nature actually being 
described by parameters in the HB/FP region. 

Toward this end, we consider Case 2, with the mSUGRA parameters set given by 

m , mi/2, A Q , tan/3, sign(fi) = 2500 GeV, 300 GeV, 0,30, +1. 

Sample sparticle masses and parameters are given in Table |3[ For these parameter choices, 
\fi\ < M2, so that the light chargino and lightest neutralino have significant higgsino com- 
ponents. The chargino mass = 113.1 GeV, and is just beyond the reach of LEP2. 
The LSP mass is = 85.6 GeV, so that the mass gap tti^ — = 27.5 GeV. The 
W\ decays via 3-body modes into Z\ff, where / and /' are SM fermions. The decays 
are dominated by the W boson exchange graphs, so that decays W\ — > Z±ff have similar 
branching fractions to W — > ff decays. 

We begin by generating e + e~ — * all SUSY particles for the signal, and generate SM 
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backgrounds using all ISAJET SM processes. We first require all events to pass the standard 
chargino pair cuts for 11 + 2j+ Et events as detailed at the beginning of Sec. 2. Next, 
following case study 4 of Ref. [p9[|, we require missing mass ijh > 240 GeV. The resulting 



signal and also background events are plotted in Fig. 12 in the E(jj) vs. m(jj) plane. 
SUSY and Higgs boson events are denoted by black dots, while SM background events 
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Figure 12: Scatter plot of SUSY signal events (black dots) and SM background (red x's) after 
standard cuts plus fa > 240 GeV cuts, in the Ejj vs. m(jj) plane. Chargino pair events occupy 
the low m(jj) region. 

are denoted by red crosses. The chargino pair events populate the cluster at low m(jj), 
since the dijet mass from chargino decay is bounded by the Wi — Z\ mass difference. If 
the chargino decays via W\ — > W Z\ — > qq'Zi, we would expect that the E(jj) distribution 
would have well-defined upper and lower endpoints that depend only on and 
(and, of course M\y), as in case study 1 of Ref. |u|]. 

In the HB/FP region, the W\ — Z\ mass gap is small, and the decay to on-shell W 
is kinematically inaccessible. We can, however, adapt this strategy by forcing "two-body 
kinematics" on these events by first selecting events in narrow bins in m(jj) and studying 



separately their E(jj) distributions. This is done in Fig. 13, where we show the E(jj) 
distribution for m(jj) bins of width 4 GeV, centered at 8, 12, 16 and 20 GeV, corresponding 
to an integrated luminosity of 100 fb _1 . This is the "data". The energy of the dijet cluster 
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is bounded by 

7 (i% - df j3 ) < E(jj) < j(E*j + (4.1) 



where = (mL + m 2 (jj) - m^^m^, p^- = - m 2 (jj), 7 = E^Jm^, (3 = 

Pjj/Ejj and £^ = \/s/2, up to energy mismeasurement errors, jet clustering, particle 
losses, bremsstrahlung and finite width bins in m(jj). The corresponding "theoretical 
predictions" shown by the smooth curve are obtained by generating a much larger sample 
of the same events and fitting this larger sample (corresponding to an integrated luminosity 
of 600 fb _1 ) to the function, 



F(E,mft.ms-,A,B,C,D)=N{l + exp 



A 



1 + cxp 



C 



Ba Emin 

"(4-2) 

where E m i n and E max are calculated for each bin in m(jj) taking the central m(jj) value 
in that bin, and input values for and m^; &E min (&E max ) ^ s the absolute value of the 
difference between E m i n (E max ) at the highest m(jj) value in each m(jj) bin and E m i n 
{Emax) for m(jj) at the center of this bin. The small contribution from the SM background 
has also been included. The parameters A, B, C and D are separately determined for 
each bin in m(jj) and serve to fit the shapes of the corresponding distributions, while N 
is adjusted to give the normalization near the maximum of Ejj distribution corresponding 
to an integrated luminosity of 100 fb . 

Next, we proceed to perform a x 2 fit to obtain and from our synthetic data 
sample, using the fitted function ( |4.2| ) for the theoretical prediction 6 for chargino and 
neutralino masses close to those for Case 2. In other words, for a grid of points in the 
(m^,m^J plane, we evaluate, 

. / F(E, m T ~r, (inp), ms (inp)) — F(E, m T -r, , m~ ) \ 
X 2 (m Wi ,m~ Zi ) = EE 1 (43) 

bins E \ yJ F ( E > m Wi> m zJ J 

where Ylbins means that we sum over all four bins in m(jj), and ^2 E denotes the sum- 
mation over all bins in Ejj and find the values of chargino and neutralino masses for 
which this quantity is minimized. These best fit values, together with the regions where 
A X 2 < 2.3 (68.3% CL) and < 4.6 (90% CL) are shown in Fig. (h| We see that it is possible 
to determine and at approximately the 10% level. 

Having determined the values of and Tn^, the next step is to examine what we 
can say about the underlying MSSM parameters fx, Mi and tan (3 that enter the chargino 
mass matrix. To determine three unknowns, we need to experimentally determine one more 
quantity which we take to be the cross section for 11 + 2j+ Et events from chargino pair 
production. Almost all the signal arises from chargino pair production if we require m{jj) < 



25 GeV and E(jj) < 100 GeV as in Fig. 12. For 100 fb -1 of integrated luminosity, we find 
1649 events, which translates to a measurement of a(e + e~ — ► W^W-j - ) = 16.5±0.4 fb (after 
all the cuts), or 2.5% statistical error. Other systematic errors will be present, although 



6 Of course, the parameters E min , UE min , E max and UE mal depend on and mg 1 via ( p~l| ) and the 

equations following that. We assume that the parameters A, B, C and D do not change. 
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Figure 13: Distribution of Ejj for £+2— jet events after standard cuts together with rfi > 240 GeV, 
with events restricted to narrow bins of rn(jj). The histograms show these distributions for the 
synthetic data sample while the solid line shows the corresponding theoretical expectation obtained 
as described in the text. 



these may be controllable by precision measurement of many SM processes. We perform a 
fit to the MSSM parameters using the values of Tn^, rn^^ and <r(W r 1 + I4 y 1 ~) as determined 
above. We scan over MSSM model parameters, using 1-loop corrected mass relations for 
and m z 1 as gi ven by ISAJET 7.69. In Fig. 15, we show the regions of a) the \x vs. M2 
plane, b) the [i vs. tan/3 plane and c) the M2 vs. tan/3 plane that are allowed at the 68.3% 
and 90% CL. In each case, we have held the parameter not shown in the plane fixed at 
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Figure 14: Fits to and and the associated error ellipses for Case 2 in the text. 

its input value. The result in frame a) clearly shows that indeed |/x| << M2, providing 
strong support that the model parameters lie in the HB/FP region, and that the LSP has 
a significant higgsino component, enhancing the neutralino pair annihilation in the early 
universe. 7 While n and M2 can be well determined (at least for this case study), it is 
also evident from the figure that a precise determination of tan/3 is not possible in this 
case. This may not be so surprising, since in the HB/FP region, SUSY scalar masses that 
depend on Yukawa couplings and hence tan/3 are so heavy that they essentially decouple 
from observable physics, and the region is relatively invariant under changes in tan (3. 

5. Summary and conclusions 

The recent constraint on the relic density of neutralinos obtained from WMAP measure- 
ments, together with earlier determinations of BF(b — > 57) and (g — 2) M select out regions 
of parameter space of the mSUGRA model. In the stau co-annihilation region, the H, A- 
annihilation funnel and in the HB/FP regions, very high values of mo and m.1/2 consistent 
with all constraints are possible: moreover, the so-called bulk region where sparticles are 
light is disfavored. These considerations motivated us to re-assess the reach of various 
collider and non-accelerator search experiments for supersymmetry. In this paper, we re- 
evaluate the reach of a -y/s = 0.5 and 1 TeV linear e + e~ collider for SUSY in the context of 
the mSUGRA model, examining for the first time the reach in the HB/FP region. We find 
that a y/s = 1 TeV LC can explore most of the stau co-annihilation region if tan p ~ 30, 
7 That \fi\ is small can presumably also be determined by studying chargino production using polarized 
beams. Note that our determination does not require this capability. 
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although along with a dilepton search, a ditau search will also be needed. The H, A- 
annihilation funnel typically extends beyond the maximum reach of a LC. In the HB/FP 
region, chargino pairs may be kinematically accessible to a LC, but the energy release in 
chargino pairs can be small, reducing detection efficiency Nonetheless, LCs should be able 
to probe much of the lower HB/FP region with standard chargino searches. In the upper 
HB/FP region, new cuts are proposed to allow signals from much of the small W\ — Z\ 
mass gap region to be observable above SM backgrounds. In this region, the reach of even 
a 500 GeV LC can exceed that of the CERN LHC! This is all the more important in that it 
occurs in a region of model parameter space which is allowed by all constraints, including 
those imposed by WMAP. 

One should also stress that the LCs reach is also 
complimentary to reach of direct dark matter search 
experiments (DDMS) even though both kinds of exper- 



parameter value (GeV) 



M 2 
Mi 
H 



m 



m 



Wi 

\v 2 



m-r 



m 



iments similarly cover much of the HB/FP region 2S]. 
The complementarity of a LC occurs for the region 
which is very close to the no REWSB border. In this 
region, the neutralino relic density is so low that DDMS 
experiments are not able to cover this part of the pa- 
rameter space even though the higgsino component of 
neutralino is significant. This region can be probed by 
experiments at a LC. 

If a supersymmetric signal is found, then the next 
obvious step will be to determine the underlying MSSM 
parameters. We have performed a case study in the low 
m i/2 part of the hyperbolic branch. In this region, we 
show that a measurement of and is possible at 
the 10% level. A measurement of the total chargino pair 
cross section to 2.5% allows a determination of MSSM 
parameters Mi and ji, although tan /3 is more difficult 
to pin down. The resulting determination of M<i and 
[i would point to a model with higgsino-like charginos 
and neutralinos. Together with absence (or low levels) 
of squark signals at the LHC, and the agreement of the 
chargino cross section with the expected s-channel con- 
tribution (pointing to heavy sneutrinos) these measure- 
ments would be indicative of an mSUGRA-type model 

in the HB/FP region. In case parameters are in the upper part of the hyperbolic branch, 
LC event characteristics may be sufficient at least to establish the production of massive 
particles, with associated decay products that are quite soft. An examination of how one 
would obtain information about the underlying scenario would be worthy of exploration. 



z 2 

m Z, 

m A 

BF{b - 



«7) 



236.5 
122.0 
121.6 
833.2 
2548.1 
2503.9 
113.1 
274.8 
85.6 
135.0 
142.2 
281.5 
2129.4 
118.8 
0.0423 
3.84 x 10~ 4 
2.3 x 10~ 10 



Table 3: Masses and parame- 
ters in GeV units for Case 2 for 
TO o, TO i/2i A( h tan/?, sigrifi = 
2500 GeV, 300 GeV, 0, 30, +1 in 
the mSUGRA model. The spectra 
is obtained using ISAJET v7.69. 
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Figure 15: Fits to fi, Mi and tan/3 from the measured values of , and cr(Wi Wj~), and 
the associated 68.3% and 90% CL regions for Case 2. The green squares denote the fitted values 
while the white squares show the corresponding input values of the parameters. 
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